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The electrostatic charging of the test mass in ASTROD I (Astrodynamical Space Test 
of Relativity using Optical Devices I) mission can affect the quality of the science data 
as a result of spurious Coulomb and Lorentz forces. To estimate the size of the resultant 
disturbances, credible predictions of charging rates and the charging noise are required. 
Using the GEANT4 software toolkit, we present a detailed Monte Carlo simulation of 
the ASTROD I test mass charging due to exposure of the spacecraft to galactic cosmic- 
ray (GCR) protons and alpha particles {'^He, "He) in the space environment. A positive 
charging rate of 33.3 e^/s at solar minimum is obtained. This figure reduces by 50% at 
solar maximum. Based on this charging rate and factoring in the contribution of minor 
cosmic-ray components, we calculate the acceleration noise and stiffness associated with 
charging. We conclude that the acceleration noise arising from Coulomb and Lorentz 
effects are well below the ASTROD I acceleration noise limit at 0.1 mHz both at solar 
minimum and maximum. The coherent Fourier components due to charging are investi- 
gated, it needs to be studied carefully in order to ensure that these do not compromise 
the quality of science data in the ASTROD I mission. 

Keywords: ASTROD I; charging; GEANT4; disturbances. 



1. Introduction 

The ASTROD I (Astrodynamical Space Test of Relativity using Optical Devices 
I) mission concept is a down-scaled version of ASTROD. The main objectives of 
ASTROD I are: to improve the precision of measurement of solar-system dynam- 
ics, solar-system constants and ephemeris; to measure the relativistic gravitational 
effects; to test the fundamental laws of space-time more precisely and to improve 
the measurement of the rate of change of the gravitational constant with time.^^ 

The basic scheme of the ASTROD I space mission is to use two-way laser interfer- 
ometric ranging and laser pulse ranging between a drag-free ASTROD I spacecraft 
in a solar orbit and deep space laser stations on Earth. The ASTROD I spacecraft 
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is 3-axis stabilized with a total mass of 300-350 kg. The mass of the payload is 
100-120 kg. The science data rate is 500 bps. The spacecraft is cylindrical with OD 
(Outer Diameter) 2.5 m and height 2 m and has its side surface covered with solar 
panels. In orbit, the cylindrical axis will be perpendicular to the orbit plane with 
the telescope pointing toward the ground laser station. The effective area to receive 
sunlight is about 5 and can generate over 500 W of power.^Hll 

The spacecraft will be launched into the solar orbit from a low earth orbit. The 
injection correction will be made using a medium-sized ion thruster. A launch on 
August 4, 2010 would provide a suitable orbit. The orbit in the X-Y plane of the 
heliocentric ecliptic coordinate system is shown in Figure 1.^2^ This solar orbit will 
initially have a period of 290 days. After two gravity-assist encounters with Venus, 
the period will be shortened to about 165 days. After about 370 days from launch, 
the spacecraft will arrive at the other side of the Sun. The spacecraft will have the 
first closest approach to Venus 107.8 days after launch with a distance of 31606.0 
km to the centre of Venus. The spacecraft crosses the Venus trajectory in front 
of Venus and gets a swing toward the Sun to achieve the Venus orbital period. 
After the first encounter, the spacecraft has the same orbit's period as Venus and 
encounters Venus again after about 1 period (224.7 days) with a closest approach 
distance of 16151.7 km from the centre of Venus.l^^^^ 
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Fig. 1. Tlie ASTROD I orbit in the X-Y plane of tlie lieliocontric ecliptic coordinate system. 

To achieve its goal, the ASTROD I residual acceleration noise target is: 

sTM) = 3 X lO--[^:^ + 30(^)^] ms-Hz-/^ (i) 

over the frequency range of 0.1 mHz < / < 100 mHz.^^Here S^J_^{f ) is the residual 
acceleration noise spectral density. It is compared to the LISA Pathfinder LISA 
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Technology Package,^ LISA'' and ASTROril noise target curves in Figure 2. 
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Fig. 2. A comparison of the target acceleration noise curves of ASTROD I, the LTP, LISA and 
ASTROD. 



The test mass is key to guarantee the drag-free condition for ASTROD I. It is 
a 1.75 kg, rectangular parallelepiped, made of an extremely low magnetic suscepti- 
bility (< 5 X 10^^) Au-Pt alloy, to minimize magnetic disturbances. The test mass 
is placed in the centre of spacecraft and surrounded by electrodes on all six sides. 
Any relative displacement of test mass to the surrounding electrodes will be capaci- 
tively detected and the spacecraft will follow it using FEEP (Field Emission Electric 
Propulsion) to ensure drag- free flight. Cosmic rays and solar energetic particles will 
easily penetrate the light shielding of the spacecraft to transfer heat, momentum 
and electrical charge to the test mass. Electrical charging is the most significant of 
these disturbances. It will result in forces on the test mass, due to Coulomb and 
Lorentz interactions, which will disturb the geodesic motion. The characteristics of 
the test mass charging process depend on the incident flux, spacecraft geometry and 
the physical processes that occur. The three main disturbances associated with this 
charge are an increase in the test mass acceleration noise, coupling between the test 
mass and the spacecraft and the appearance of coherent Fourier components in the 
measurement bandwidth. To limit the acceleration noise associated with Coulomb 
and Lorentz forces to meet the ASTROD I noise requirement, the test mass must be 
discharged in orbit. Our previous work predicted the charging rates for ASTROD I 
test mass from galactic cosmic rays at solar minimum using a simplified geometry, 
and using these predictions, estimated the magnitude of disturbances associated 
with charging.l^ In this paper, we present the detailed calculation of the ASTROD 
I net test mass charging rate and shot noise, due to cosmic rays at solar minimum 
and solar maximum, with a more realistic geometry model. Based on these results, 
we estimate the magnitude of acceleration noise, stiffness and the coherent signals 
associated with charging. 
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2. Modelling the Charging Process 
2.1. Radiation environment model 

We have simulated the fluxes of the 3 most abundant primaries, proton, '^He and ''He 
primary particles which represent approximately 98% of the total cosmic ray flux. 
Near-Earth cosmic ray spectra were adopted, as used in similar LISA simulations. 
These are shown in Figure 3.^^ The primary particles are generated from points 




GCR Energy (GeV) 



Fig. 3. Differential energy spectra for cosmic ray protons and He nuclei at near earth orbit. For 
each species, the upper curve indicates the solar minimum spectrum, the lower curve indicates the 
solar maximum spectrum. 

sampled uniformly from a spherical surface of 3000 mm diameter, encompassing 
the whole ASTROD I geometry model. 

The total time T for bombardment by Nq primaries is given by T = A^o /F ■ ttB? , 
where F is the integral flux (unit: particles/cm^ /s) for each species at solar minimum 
or solar max imu m, and R = 1500 mm. The A^*'' particle event occurs at time 
t = N ■ T/Nq^^^ The effects on charging of other particle species (C, N, O, e~) are 
determined separately, based on a LISA study.^^ Work is under way to study the 
impact on the charging disturbances from variations in the incident flux due to the 
ASTROD I heliocentric position changes. 

2.2. Physics model 

Test mass charging depends heavily on the physics processes that occur during the 
passage of the particles through matter and the geometry model used in simula- 
tion. The GEANT4 toolkit employs Monte Carlo particle ray-tracing techniques to 
follow all primary and secondary particles. Due to their high energy and hadronic 
nature, cosmic rays can produce the complex nuclear reactions which have large 
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final-state multiplicities, producing many secondary particles. A low energy thresh- 
old of 250 eV was imposed for secondary particle production in our simulation. The 
physics processes simulated include electromagnetic, hadronic and photonuclear in- 
teractions. Fluorescence and non-radiative (Auger) atomic deexcitation have been 
implemented. The hadronic physics is mainly implemented by elastic and inelastic 
scattering processes. The inelastic reactions were based on the LEP (Low Energy 
Particles) and HEP (High Energy Particles) parameterized models. The inelastic re- 
actions also use evaporation models to treat the deexcitation of nuclei with A > 16, 
comprising gamma emission, fragment evaporation (p, n, a, ^H and '^H) and fission 
of heavier residual nuclei. A variety of decay, capture and annihilation processes has 
also been included in our physics processes list based on LISA GEANT4 model of 
Araujo et al.l^ The charging potential of several additional physics processes, such 
as the kinetic emission of very low energy electrons which has no t b een modeled in 
the present simulation, has been assessed based on LISA studies.^^^ 

2.3. Geometry model 

The basic payload configuration of ASTROD I is sketched in Ref. 1. The geometry 
model built for ASTROD I using GEANT4 in present work is sketched in Figure 
4 and Figure 5. Table 1 lists the dimensions and weight of main constituent parts 
of the ASTROD I geometry model. Table 2 lists the composition and density of 
materials used in the model. 

The cylindrical structure of the spacecraft consists of a layer with diameter 
2.5 m, height 2 m and thickness 10 mm made of CFRP (Carbon Fibre Reinforced 
Plastic) honeycomb. The top and bottom of the spacecraft are covered by the upper 
deck and lower deck. To prevent sunlight from striking the inside directly and to 
reduce the temperature perturbation inside, all surfaces of spacecraft are covered 
by a thermal shield consistent of five layers of materials (face sheet, honey comb 
core, face sheet, foam, face sheet). The edges of upper deck and thermal shield are 
shown as large ellipses in Figure 4. The inner lower deck is shown as the grey part in 
Figure 4 and Figure 5. The payload structure is used for shielding the optical bench, 
inertial sensor and primary telescope etc. The primary telescope , wh ich collects the 
incoming light is a 500 mm diameter f/1 Cassegrain telescope.^^ Some 30 boxes 
represent the components above ~ 0.1 kg based on the payload configuration.^I] The 
mass of the spacecraft and payload are estimated as 341 kg and 109 kg in this study. 
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Fig. 5. The overhead view of the GEANT4 geometry model for ASTROD I. The six black boxes 
are laser heads and the other white boxes represent the components above ~ 0.1 kg. 
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Table 1. Dimension and weight of the main constituent parts of the ASTROD I geometry model 



Constituent part 


Material 


Weight 
(kg) 


Dimension (mm) 


Comment 


Solar Panels 


Scell 


~61.5 


Tube: radius 1250; 
height 2000; thickness 0.5 


covers the side 
surface of the 
spacecraft 


Spacecraft 


CFRP 


~7.8 


Tube: radius 1250; 
height 2000; thickness 10 


contains the 
payloads 


Thermal Shield 


CFRP 

Al-Honeycomb 
foam 


~100.3 


Tube: radius 1240; 
height 2000; 
thickness 41.8 


covers the 
spacecraft 


Upper Deck 


CFRP 

Al-Honeycomb 


~15 


Cylinder: radius 1240; 
thickness 30 


top of the 
spacecraft 


Lower Deck 


CFRP 

Al-Honeycomb 


■^15 


Cylinder: radius 1240; 
thickness 30 


above the 
MLI-blanket 


TV (TT T 1- 1 1 J. 

MLi- blanket 


MLImat 


~6.8 


Cylinder: radius 1240; 
thickness 1 


bottom of the 
spacecraft 


Payload Shield 


CFRP 


~15.6 


Tube: radius 300; 
length 1000; thickness 5 


Shields the 

optical bench, 
incrtial sensor 


Optical Bench 


ULEglass 


~5 


Rectangular: length 350; 
width 200; height 40 


contains the 
Ti-house 


Ti House 


Ti Alloy 


~1.9 


lube: radius 62.5; 
height 224; thickness 5 


located in the 
optical bench 


Mo House 


Molybdenum 


~2.5 


Cube: length 75 


contains the 
test mass 


Test Mass 


AuPt Alloy 


~1.75 


T~> A. 1 1 j_i rri 

Rectangular: length 50; 
width 50; height 35 




Primary 
Telescope 




'•^6.4 


Dish: radius 250; 
thickness 10.5 


collects the 
incoming light 


Shield/Mouutius 
1'la.le 


CFHP 


~1.6 


Cyliiidor: radius 175: 
t.liickiiess 10 




Telescope Shield 


CFRP 


~12.6 


Tube: radius 285; 
length 860; thickness 5 


Shields the 
secondary mirror 


PCDU 


A16061 


~15.9 


Rectangular: length 350; 
width 200; height 300 


Power Condition- 
ing and 

Distribution Unit 


Transponder 


A16061 


~3.5 


Rectangular: length 220; 
width 184; height 178 




CPS 


A16061 


~15.9 


Rectangular: length 240; 
width 356; height 140 


Centralised 
Processor System 


Interferometer 
Electronic Boxes 


A16061 


~3.5 


Rectangular: length 200; 
width 200; height 150 




Gyroscope 


A16061 


~1 


Cylinder: radius 42.5; 
height 89 




RFDU 


A16061 


~1 


Rectangular: length 160; 
width 60; height 80 


Radio Frequency 
Distribution Unit 
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Table 2. The composition and density of materials used in the ASTROD I geometry 
model (CFRP;Carbon Fibre Reinforced Plastic) 



Material 


Composition (by weight) 


Density (g/cm'^) 


Vacuum 


gas 


1.0x10^^^ 


A16061 


Al(98%), Mg(l%), Si(0.6%), Fe(0.4%) 


2.70 


Al Honeycomb 


Al(98%), Mg(l%), Si(0.6%), Fe(0.4%) 


0.05 


MLImat 


H(4.1958%), C(62.5017%), 0(33.3025%) 


1.40 


Ti Alloy 


Ti(90%), Al(6%), V(4%) 


4.43 


AuPt Alloy 


Au(70%), Pt(30%) 


19.92 


ULE Glass 


SiGlass(92.5%), TiGlass(7.5%) 


2.21 


foam 


C(90%), H(10%) 


0.05 


SiGlass 


O: 2, Si: 1 


2.20 


TiGlass 


O: 2, Ti: 1 


4.25 


SHAPAL 


Al: 2, N: 1 


2.90 


SiC 


Si: 1, C: 1 


3.10 


Scell 


Si 


7.82 


Molybdenum 


Mo 


10.22 


Gold 


Au 


19.32 


CFRP 


C 


1.66 


carbon 


C 


2.10 


CFRP Honeycomb 


C 


0.05 



A 50 X 50 X 35 mm^ test mass is at the centre of the spacecraft. The test 
mass is housed inside capacitance sensors located in optical bench mounted behind 
the telescope. The test mass is surrounded by sensing and actuation electrodes 
lodged in a molybdenum housing. A 0.3 [xm gold layer is plated on the entire inner 
surface of the sensor housing. The assembly is accommodated in a titanium vacuum 
(< 10 uPa) enclosure. The gap between test mass and electrodes along X axis or Y 
axis is 4 mm; that along Z axis is 2 mm.^l^ The GEANT4 model for inertial sensor 
is shown in Figure 6. 




Fig. 6. ASTROD I inertial sensor model implemented in Geant4. The test mass, located at the 
centre of the figure, is surrounded by sensing electrodes (white) and injection electrodes (grey). 
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3. Charging Results 

We have run 6 independent GEANT4 simulations to determine the charging of 
the ASTROD I test mass by cosmic ray protons, ^He and ^He, at solar minimum 
and maximum. In total, about 8,500,000 events were simulated. The details of each 
event that resulted in test mass charging were recorded, including the event time, 
net charge deposited on the test mass and the energy of the primary. 

3.1. Charging simulation results 

The variation of the net test mass charge with time, due to GCR proton, ^He 
and *H{^ fluxes are shown respectively in Figure 7, Figure 8 and Figure 9 at solar 
minimum and maximum. The straight lines in these figures correspond to least 
squares fits of the simulated data, giving mean net charging rates attributable to 
the proton, ^He and "'He fluxes. The proton flux is responsible for positive charging 
rates of 26.5 ± 0.5 e+/s at solar minimum and 9.0 ± 0.5 e^/s at solar maximum; 
The "^He flux is responsible for positive charging rates of 0.8 ± 0.05 e^/s at solar 
minimum and 0.3 ± 0.05 e+/s at solar maximum; The ^He flux is responsible for 
positive charging rates of 6.0 ± 0.20 e+/s at solar minimum and 2.4 ± 0.20 e+/s at 
solar maximum. The uncertainties quoted are only associated with the Monte Carlo 
fluctuations. Our simulation indicates that ~ 97% of the charge accumulated comes 
from primary cosmic ray protons and *He at both solar minimum and maximum 
and all three fluxes lead to positive charging of the test mass. The proton flux 
dominates these rates. However, ''He, which constitutes only 8% of the total cosmic 
rays flux, is responsible for ~ 18% of the test mass charging at solar minimum and 
20% at solar maximum. 
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Fig. 7. Tile cliarging timeline for protons at solar minimum and maximum. The straight line is 
a least squares fit. 
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Fig. 9. The charging timeline for ^He at solar minimum and maximum. The straight line is a 
least squares fit. 



Two histograms of the net charge deposited in an event are given in Figure 
10 and Figure 11, for the proton data set, showing that most events result in the 
transfer of one unit of charge. The effects of the positive and negative chargings 
cancc4 to some extent. An imbalance in these currents gives rise to the net positive 
charging rate. 



I 



-20 -15 -10 -5 5 10 15 20 

Net Charge Deposited in Envent 



Fig. 10. Histogram of the net charge deposited in an event, for incident protons at solar minimum. 
The total number of proton events simulated was 2,290,000. 
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Fig. 11. Histogram of the net charge deposited in an event, for incident protons at solar maximum. 
The total number of proton events simulated was 4,000,000. 



The charging rate is plotted as a function of primary energy in Figure 12. The 
low energy cut-off is due to the shielding provided by the spacecraft, which prevents 
incident protons with energies below ~ 100 MeV from charging the test mass. At 
solar minimum, the most significant charging mechanism is primary cosmic ray 
particles stopping in the test mass. This occurs mainly for protons of energy between 
~ 100 - 720 MeV. Protons with energies in excess of ~ 720 MeV have sufficient 
energy to traverse the distance through the spacecraft to the test mass and the 
longest path through the test mass, without being stopped. This explains the peak 
observed in this energy interval in Figure 12 at solar minimum. The scenario at 
solar maximum is distinct: a peak is visible at higher energies because the primary 
proton flux peak shifts towards higher energy. 
3.5 p 

3.0 - 

2.5 - 

2.0 - 

a 

I 1.5 - 
3 1.0 - 
i 0.5 - 
I 0.0- 

-0.5 - 

-1.0 - 

GCR Energy, MeV 

Fig. 12. Charging rate of test mass as a function of primary proton energy at solar minimum and 
at solar maximum. 
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In addition to the Monte Carlo uncertainty, we added an error of ±30% in the 
net charging rates to account for uncertainties in the GCR spec tra, ph ysics models 
and geometry implementation. Further, based on LISA studieJim, a potential 
contribution to the charging rate of 28.4 e"'"/s at solar minimum and of 17.0 e^/s 
at solar maximum from kinetic low energy secondary electron emission should be 
considered; the effect of cosmic ray fluxes of particle species not included in this 
simulation is expected to incre ase the net charging rate by ~ 4.2% at solar minimum 
and 7.3% at solar maximu nP. 



3.2. Charging noise 

Following to Ref. 12, the charging flux is considered to be made up of independent 
currents, Iq, each composed solely of charges qe {q=+l for protons), with shot noise 
of single-sided spectral density Sq — -^2ge/g, where e is the magnitude of electron 
charge. The total noise, Sr, is then given by the quadrature sum of Sq, over all values 
of q. Considering the Monte Carlo currents alone, Sq = 17.6 es~^Hz~^/^ at solar 
minimum and Sq = 9.6 es~^Hz~^/^ at solar maximum. Based on the LISA studjC^, 
low energy secondary electron emission is estimated to contribute an extra 10.3 
es~^Hz~^/^ at solar minimum and 8.0 es~^Hz~^/^ at solar maximum; the charging 
noises from other species (C, N, O, e~) are estimated to be 8.9 es^^Hz^^/^ at solar 
minimum and 3.7 es^^Hz^^/^ at solar maximum. Integrating in the time domain 
gives the charging fluctuations at frequency /, 

SqU) = SR/2nf. (2) 

The charging rate and noise contributions from the different sources mentioned 
above are summarized in Table 3. Given the results in Table 3, by adding the con- 
tributions of all the independent sources, we estimate the total worst case charging 
rate to be 73.8 e^/s at solar minimum and 33.8 e+/s at solar maximum. The total 
worst case noise is 22.3 es~^Hz~^/^ at solar minimum and 13.0 es^^Hz^^/^ at solar 
maximum. 



Table 3. The charging rate and noise contributions from different sources. 



Source 


charging 


rate{e+/s) 


charging 


noise(e/s/Hz-'-''^) 




min 


max 


min 


max 


P 


26.5 


9.0 


15.9 


8.6 


3He 


0.8 


0.3 


2.5 


1.6 


"He 


6.0 


2.4 


7.2 


3.9 


Secondary Electron 


28.4 


17.0 


10.3 


8.0 


Other Species{C, N, O, e") 


1.4 


0.9 


8.9 


3.7 


Uncertainty 


10.0 


3.5 
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4. Acceleration Disturbances 
4.1. Coulomb noise and stiffness 

The charge-dependent Coulomb acceleration a^j^ in direction k is given by: 



g2 dCr QVt OCt q ^,,da 



2mC2 dk mCr dk mCr ^ ' dk ' ^ ' 

^ i—l 

The first two terms in equation (3) are dependent on the overall sensor geometric 
symmetry, through dC'T/dk, and the third term is dependent on the symmetry of 
the sensor voltage distribution. The corresponding acceleration noise, Saqk^ due to 
random fluctuations of the test mass position relative to the spacecraft, 6k, of the 
potentials of the conductors that surround the test mass, SVi, and of the test mass 
free charge 6Q, is given by: 

z— 1 ^ 

where k is a displacement in direction fc; m is the mass of the test mass; Q is the free 
charge accumulated on the test mass; C'ij is the capacitance between conductors 
i and j which surround the test mass; Vt is the potential to which conductor i is 
raised; Ct = X^ila^ Ci,Af is the coefficient of capacitance of the test mass, which 
is defined as the TV*'' conductor, with potential Vat = Q/Ct + Vt, and Vr = 
^Y^C,. V."^ The estimates for acceleration noise have assumed typical 
parameter values for the ASTROD I mission: Q was taken as the amount of charge 
accumulated in 1 day, assuming a net test mass charging rate of 73.8 e~''/s at solar 
minimum and 33.8 e+/s at solar maximum, which corresponds to the Monte Carlo 
rate, with error margins, estimated contributions from particle species not included 
in the Monte Carlo model and the potential contribution from kinetic low energy 
secondary electron emission, that is likely to almost cancel in the actual sensor,^!^ 
added; m =1.75 kg; mean voltages on opposing conductors Vi = 0.5 V; the potential 
difference between conductors on opposing faces of the sensor compensated to 10 
mV; the asymmetry in gap across opposite sides of test mass is 10 /J,m; capacitances 
and capacitance gradients were calculated using parallel plate approximations: Ct 
= 53 pF; Vt = 0.5 V; position noise Sk = 1x10^^ mHz~^/^; voltage noise SVt = 
1x10-* VHz-i/2 and charge noise SQ = 4.6x10 CHz which includes, as for 
the charging rate, the unmodelled contributions. These Coulomb acceleration noise 
due to the test mass charging are listed in Table 4. The total noise estimated here 
is a factor of ~ 30 less than the ASTROD I acceleration noise target. The stiffness 
associated with test mass charging, Sqk, is given by Sgk = —m ■ dagk/dk. These 
acceler atio n noise figures are lower than the results liberally estimated by Shiomi 
and nPI because of the smaller charging rate and charging noise we obtained 
here. The requirements on Coulomb noise, stiffness and other associated noises for 
ASTROD I in Ref. 14 are satisfied. 
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Table 4. The magnitude of charging disturbances for ASTROD I at solar minimum and 
maximum, at frequency = 0.1 mHz. 



Solar activity 




Coulomb 


noise 




Lorentz noise 


Stiffness 




(xlO -"-^ms 






(xlO-iSj^s-^Hz-O-^) 


(xl0-**s-2) 




Sk 




5Q 


total 






minimum 


1.06 


2.18 


1.32 


2.80 


2.80 


-1.52 


maximum 


0.49 


1.08 


0.77 


1.40 


1.30 


-0.69 



dk: displacement noise; SVi: voltage noise; 6Q: charge noise. 



4.2. Lorentz noise 

Lorentz effects arise from the motion of the test mass through the interplanetary 
magnetic field, and its residual motion through the field generated within the 
spacecraft, . The test mass will be housed in a conducting enclosure, which will 
reduce the effect of the interplanetary field, via the Hall effect, with efficiency t]. 
Hence, to first order, the Lorentz acceleration noise, a^^, is given by: 

m'K)' = ivQV.SB^)^ + {iiQ5V,B,f + {Q5V,B^f + {r,6QV,Bj^. (5) 

where Vj is the speed of the test mass through the interplanetary field; SVj and 
SVs are the magnitudes of random fiuctuations in the test mass velocity through 
the interplanetary field and relative to the spacecraft, respectively and SBj gives 
the magnitude of fiuctuations in the interplanetary field. ' a^^ also increases with 
decreasing frequency, and is estimated to be ^ 2.8x10"^^ ms~^Hz^^/^ (0.1 mHz) 
at solar minimum and ~ 1.3x^^^ nis^^Hz^^/^ (0.1 mHz) at solar maximum, which 
is a factor of ~ 30 below the ASTROD I acceleration noise target. We have assumed 
that Q — 73.8 e+/s at solar minimum and 33.8 e+/s at solar maximum, as in section 
4.1; rj = 0.1; \7 = 4x10"^ m/s; SV,^ 4.78x10-1^ ms-^Rz-^/^; 6V^ = 6.28x10"" 
ms-iHz-i/2; = 9.6x10-6 T; \ 6B^ \ = 1x10"^ THz^i/^. 1.2x10^^ T 

(this is a conservative estimate of the field at 0.5 AU, used to give the worst-case 
noise, for the ASTROD I orbit) and \ 6B^ \ ^ 1.2x10"^ THz^^/^ xhe estimates of 
acceleration noises and stiffness due to the Coulomb and Lorentz effects are listed 
in Table 4. 

4.3. Coherent Fourier components 

The charging of the test mass may also result in the appearance of unwanted, 
coherent Fourier components in the ASTROD I measurement bandwidth through 
Coulomb and Lorentz interactions, due to the time dependence of the amount of 
charge accumulated on the test mass. It is shown that the signals associated with 
Lorentz interactions are expected to fall below the ASTROD I test mass residual 
acceleration noise, but the signals from Coulomb interactions may exceed both the 
instrumental noise over a fraction of the bandwidth, and may not be eliminated by 
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daily discharging of the test mass. The free charge on the test mass at time t can 
be expressed as follows: 

Q{t) « '^t. (6) 

where t is the time for which the test mass has been allowed to charge and Q is the 
mean charging rate. Here, Q is assumed to be constant. Substituting Q{t) « Qt into 
the expressions for the Coulomb and the Lorentz accelerations gives the termJlSl 

T 

e.(t).e...-^|, (8) 

lk{t) = <i>kt='^{V, X B,)-k, (9) 

where /fc(i) is dependent on the overall sensor geometric symmetry; efe(i) is depen- 
dent on the symmetry of the sensor voltage distribution; lk{t) is caused by the inter- 
planetary magnetic field. The Fourier transforms of these signals with fixed-interval 
discharge will be described by a series of sine functions {sinc{x) = simrx/Trx). The 
equivalent one-sided power spectral density for these coherent signals is given by 
Ppj, = 2FT'^ /t, where FT = Fourier transform of the signal and r is the length of 
the data sample time. 

Implementing the parameter values given in section 4.1 and 4.2, taking the mean 
charging rate as constant and assuming that the test mass is discharged once every 
24 hours (as described in Ref. 15), the spectral densities of efe(t) and lk{t) are 

estimated. The curves in figure 13 and figure 14 trace Pp^i^f), PskU) ^ife(/)' 
at the primary peaks of the sine functions, where Fk{f), Ek{f) and Lk{f) are the 
Fourier transform of fk{t), ek{t) and lk{t). fk{t) and ek{t) exceed the ASTROD I 
acceleration noise limit of 3 x 10"" [0.3 mHz// -I- 30(//3 mHz)^] ms^^Hz-i/^ in a 
fraction of frequency bins in the frequency range of 0.1 mHz < / < 4 mHz at solar 
minimum (see Figure 13); only ek{t) exceeds the limit in a fraction of frequency bins 
in the frequency range of 0.1 mHz < / < 2 mHz at solar maximum (see Figure 14). 
These effects are more severe as frequency decreases. Several schemes could be used 
to minimize a potential loss of the ASTROD I science data, including continuously 
discharging the test mas s, m inimizing sensor voltage and geometrical offsets and 
through spectral analysis!--' 

Variations in, for example, the mean charging rate, could result in these signals 
exceeding the noise target in a larger fraction of the bandwidth. Hence, variations 
in these signals need to be studied carefully as they will influence the accuracy with 
which the solar-system and relativistic parameters can be determined. 

A charge management system has been developed by Imperial College London, 
for LISA Pathfinder. This system has been extensively tested, both via simulations 
and laboratory tests. A similar system could easily be used for ASTROD I. 
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Fig. 13. The curves trace the spectral densities at the primary peaks of the sine functions at 
solar miniinuin, of the coherent Fourier components, for t=1 year: ej.(t) is given by the dotted 
Une, /fj it) is given by the dashed line and /j. (t) is given by the dashed-dotted line. The bold full 
line gives the ASTROD I acceleration noise limit. 
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Fig. 14. The curves trace the spectral densities at the primary peaks of the sine functions at 
solar maximum, of the coherent Fourier components, for t=1 year: ef^{t) is given by the dotted 
line, fk{t) is given by the dashed line and lk{t) is given by the dashed-dotted line. The bold full 
line gives the ASTROD I acceleration noise limit. 



5. Conclusion 

The charging of the ASTROD I test mass by cosmic ray protons and alpha particles 
(^He and ^He) has been simulated using the GEANT4 toolkit at solar minimum 
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and maximum. The Monte Carlo model predicted a net charging rate of ~ 11.7 
e"'"/s at solar maximum, rising to 33.3 e^/s at solar minimum. Although the proton 
flux is the dominant charging flux, *He, which constitutes only 8% of the total 
cosmic ray flux, is responsible for ~ 18% and ~ 20% of this rate at solar minimum 
and maximum, respectively. We have also included an additional net charging rate 
contribution due to particle species that were not included in the Monte Carlo 
model, and a potential charging mechan ism, d ue to kinetic low energy secondary 
electron emission based on LISA studies.HSEU There is an additional uncertainty 
of ± 30% in the net charging rate, due to uncertainties in the cosmic ray spectra, 
physics models and geometry implementation. A recent, preliminary, comparison of 
a simplified GEANT4 simulation and actual GP-B charging rates, indicate 45% 
agreement (GP-B simulation is la rger than experimental value by 45%), reinforcing 
confidence in these predictions.^^^EH 

The ASTROD I acceleration noise limit is 10"^'^ ms^^Hz^^/^ at 0.1 mHz, which 
is less stringent than the LISA requirement. The magnitudes of the Coulomb and 
Lorentz acceleration noise associated with test mass charging increase with decreas- 
ing frequency. At the lowest frequency in the ASTROD I bandwidth, 0.1 mHz, the 
estimates of the Coulomb and Lorentz acceleration noise are both well below the 
acceleration noise target. These results agree, to within 30%, with those from our 
earlier studj^, which was based on a simple geometry model. The variations in the 
test mass charging rate will alter the spectral description of the coherent Fourier 
components. Hence, further work is needed to ensure that these do not compromise 
the quality of the science data of the ASTROD I mission. 

The charging process of the ASTROD I test mass by SEPs (Solar Energetic 
Particles) has also been simulated using the GEANT4 toolkit, and the charging 
rate is much larger than the values due to GCR (Glactic Cosmic Ray) proton flux 
at solar maximum and solar minimum !2Q|-[22| jjowever, the charge management 
hardware described in reference 18 could be used to discharge a test mass even 
during solar events, provided safe operation could be ensured. 

The effect of cosmic ray fluxes of particle species not included in the Monte 
Carlo simulation needs to be verified for the ASTROD I geometry. According to 
the recent work by Grimani et al.,'^ the charging rate in each LISA test mass 
induced by primary and interplanetary electrons is comparable (absolute value) to 
that released by the nuclei of the C, N, O group at solar minimum. The acceleration 
disturbances due to charging of the ASTROD I test mass by electrons are also under 
study. Further, we will evaluate the variation in the ASTROD I test mass charging 
rate over the orbit over the solar cycle, including a detailed study of SEP events, 
and its variation due to modulation of cosmic ray flux over the ASTROD I orbit. 
For this, a SCoRE (Solar And Cosmic Ray Physics And The Space Environment: 
Studies For And With LISA) study along the line of Shaul et al.— would be useful. 
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